Ecology Letters

The shape of density-dependence governs species persistence, and ecosystem stability. Yet, whether per-capita growth declines sublinearily, or superlinearily with density remains hotly debated. Growth rates across the tree of life have been shown to decline sublinearly with density, whereas theory founded on resource competition predicts the opposite. Here, we resolve this discrepancy and show that sublinearity can readily emerge from geometric constraints on consumer interactions. By linking inter individual spacing, movement and interference rates, we derive two limiting-interference regimes, one of which the well-mixed limit recovers the form of classic Beddington DeAngelis interference response. We then developed an individual-based model from first principles which reproduces the derived sublinearity response, and further use empirical data from published consumer-resource experiments that also bears the signature of sublinear density-dependence. Further, embedding the interference mechanisms underlying the emergence of sublinear density-dependence in coexistence theory opens a new regime for species coexistence where classical theory fails to predict. Our framework indicates that non-consumptive interactions are not merely a correction to resource competition but might be a distinct axis along which diverse communities may potentially coexist.

Most evidence of the consequences of climate change in natural systems has focussed on shifts in mean temperature (1,2), but the effects of extreme climatic events (ECEs) remain far less understood. This is particularly true for very severe ECEs that may occur only once every few decades. Understanding the consequences of these severe events for natural populations is nonetheless critical, since their frequency is predicted to rise under current climate change (3). Here we combine a unique long-term dataset spanning almost five decades of breeding (>20,000 events) and morphological data (>120,000 observations) in adult and nestling great tits (Parus major) and blue tits (Cyanistes caeruleus) with fine-scale temperature records to examine the effects of an unprecedented heatwave in May 2026 on breeding success and morphology. Average temperature during the heatwave (22-29 May 2026) was 7.85 C above the historical record, reaching +10.5 C (+4.32 SD) at its peak (25-26 May). These record-breaking temperatures significantly reduced adult breeding success and nestling bmass relative to expectation in the absence of a heat-wave. Given the heatwave was widespread (Fig. 1A), our findings from a single, exceptionally well-studied population are likely to generalise to other species exposed to the same event, providing key evidence that severe ECEs can substantially harm wild populations.

Ecologists remotely track movement steps of animals (e.g., via global positioning systems) and use step selection functions to study the effect of environmental factors upon their movement decisions. Constructing such functions requires pairing each observed step with some number of unobserved but feasible comparison steps. Larger numbers of comparison steps generally yield better estimates but also incur potentially challenging computational demands. Thus, it is important to determine an appropriate number of comparison steps. No established guidance exists for this decision. Here, we use simulated tracks to assess how many comparison steps are needed, fitting each set of steps to a conditional logistic regression model. We monitor errors in estimated effects for several classes of tracks, identifying the number of comparison steps for which mean relative absolute error in estimated effects is consistently low. By this criterion, 32 comparison steps per observed step are needed for our primary class of simulated tracks. Tracks in more homogeneous landscapes, tracks with shorter mean step lengths, or shorter tracks generally require more comparison steps (ranging from 64 to 128 per observed step) to achieve the same level of accuracy. Longer tracks generally require fewer comparison steps (16 per observed step). These results clearly demonstrate that the number of comparison steps influences how well step selection functions estimate covariate effects and provides initial direction in a research area that currently lacks quantitative guidance. Movement ecologists should take care when selecting the number of comparison steps paired with each observed step because those decisions matter.

Rapid evolution allows populations to persist in environments where they would otherwise go extinct. This phenomenon, known as evolutionary rescue, is typically studied in the framework of biological evolution, yet adaptive traits can also arise and spread through cultural evolution. The present study developed a stochastic eco-evolutionary model to compare rescue probabilities through biological and cultural evolution. Transmission bias governed the rescue probability under cultural evolution by setting how readily a rare adaptive trait was copied. Conformity bias suppressed population persistence because a rare trait was the least likely to be copied. Content bias toward the adaptive trait enabled evolutionary rescue when social learning was rapid, but it typically yielded a lower rescue probability than biological evolution. Only anticonformity bias, together with a high social learning rate, exceeded the rescue probability of biological evolution by enabling the adaptive trait to be established more rapidly. These results demonstrate that transmission bias alters the demographic consequences of cultural evolution and highlight the importance of transmission processes in evolutionary rescue theory. Understanding how adaptive behaviours are socially transmitted may also improve predictions of animal population persistence and inform conservation efforts in rapidly changing environments.

Squid abundance has been reported to increase globally between 1970 and 2010. This increase has been hypothesized to result from two primary factors: the loss of top predators due to overfishing and rising ocean temperatures. The decline in apex predators may lead to the expansion of squid populations either through reduced predation pressure or diminished competition with juvenile predators. Concurrently, increased temperatures could enhance the somatic growth rates of squid, thereby accelerating their population growth. However, empirically disentangling the impacts of predator loss and temperature on squid biomass remains challenging, especially in a food-web context. In this study, we used a size- and trait-based model of upper trophic levels that resolves the ecosystem structure -- biomass and trophic interactions of fish and squid -- for varying depth, temperature, and secondary production, to investigate two hypotheses of the historical expansion of squid, i.e., the effects of predator depletion from fishing and rising temperatures on squid biomass. Our model reveals that intensified fishing of squid predators -- specifically large demersal fish in shelf systems and large pelagic fish in open oceans -- leads to a slight increase in squid biomass. Conversely, elevated temperatures are associated with a decline in squid biomass. This temperature-driven reduction in biomass is attributed to an increased metabolism of squids beyond the available food supply. If historic overfishing on large marine predators continues to be curtailed, we expect a corresponding reduction in global squid biomass and fisheries potential, which could be further exacerbated by rising temperatures.

Urbanization is reshaping freshwater ecosystems, with well-documented effects across gradients of land-use change, hydrologic alteration, and habitat degradation. However, how biodiversity is organized among neighboring urban aquatic habitats that differ in hydrologic connectivity, disturbance transmission, residence time, management history, and opportunities for species movement is often less clear. This creates a challenge for interpreting urban fish communities at local scales as species occurrence may reflect both contemporary habitat filtering and historical contingencies including native persistence, interbasin transfer, stocking, and nonindigenous introductions. Here we use eDNA detections, historical records, phylogenetic information, and species trait data to investigate the fish assemblages of the Charlotte metropolitan region. We detect a highly mixed fauna that also depicts a strong signature of structured biodiversity profiles across taxonomic, phylogenetic, functional, and life-history dimensions between habitat types. In particular, bounded habitats contained assemblages with larger-bodied species that are fecund and faster to reproduce relative to free-flowing habitats. Species-level occurrence models did not support a simple trait-by-habitat rule. Instead our results demonstrate that urban aquatic habitats can sort historically mixed regional species pools into predictable assemblage-level life-history profiles while simultaneously retaining signatures of evolutionary and historical biogeographic contingency.

Viral infection is a major yet poorly quantified driver of microbial interactions and biogeochemical fluxes in the ocean. In diatoms, which are key contributors to marine primary production, the extent to which viruses reprogram host cell metabolism and alter elemental cycling remains largely unresolved. Here, we investigated how infection by a lytic single-stranded RNA virus reshapes carbon (C) and nitrogen (N) fluxes in the ecologically relevant nanoplanktonic diatom Mediolabrus comicus. Using controlled infection experiments coupled with flow cytometry, electron microscopy, PAM fluorimetry, and stable isotope probing, we resolved infection-driven changes from the population to the cellular scale. Infection induced rapid optical shifts and cellular reorganization, including the formation of membrane-bound viral replication compartments. These changes coincided with early impairment of plastidial functions, as shown by disruption of photosystem II functionality and a concomitant decline in photosynthetic carbon fixation. In contrast, nitrogen uptake was maintained and strongly enhanced during late stages of infection, indicating sustained resource acquisition to support viral replication. This decoupling led to dynamic changes in cellular stoichiometry and, overall, to substantial reductions in population-level carbon and nitrogen assimilation due to growth inhibition. Together, these findings demonstrate that diatom RNA virus infection reshapes host carbon and nitrogen metabolism, with cascading effects on elemental cycling. Our results identify diatom RNA viruses as important drivers of marine biogeochemical processes, with implications for primary production and the fate of organic matter in the ocean.

Anterograde and retrograde signalling establish bidirectional communication between the nucleus and chloroplasts. Retrograde signals from chloroplasts regulate nuclear gene expression while anterograde signals from the nucleus coordinate chloroplast development and maintain cellular homeostasis. How this bidirectional signalling framework extends beyond locus-specific regulation to shape the global spatial organization of nuclear chromatin across tissues remains unclear. Although anaesthesia can alter chromatin organisation, the role of chloroplast dysfunction in these changes remains unclear. Here, we investigate how chloroplast dysfunction and anaesthesia influence euchromatin and heterochromatin organisation in Solanum lycopersicum seedlings across tissues with contrasting photosynthetic competence. Using confocal and super-resolution radial fluctuation (SRRF) imaging with quantitative multiparameter analysis, we identify distinct, tissue-specific chromatin responses to chloroplast disruption and anaesthesia. Notably, anaesthesia induces distinct spatial chromatin changes across tissues that are independent of chloroplast dysfunction, suggesting a direct nuclear response to anaesthesia rather than a chloroplast-mediated retrograde effect. These findings highlight chromatin topology as a potential quantitative biomarker of cellular disruption and provide a framework for investigating anterograde chloroplast-nucleus coordination and stress-responsive nuclear organisation in plants.

1. Forest canopies harbour strong microclimatic gradients that shape plant performance, species interactions and ecosystem processes. Yet, despite renewed interest sparked by global change, forest canopies remain difficult-to-access experimental spaces. 2. With the goal to expand access to tree canopies as experimental arenas, we designed, built, and tested TreeTOP, a standardized experimental platform that opens canopy space for manipulative ecological experiments, specifically with potted plants. TreeTOP features lightweight aluminum frames placed in mature tree canopies non-invasively, allowing potted plants to be placed in three different heights, ground level, shade canopy, and sun canopy. 3. We implemented TreeTOP using two contrasting infrastructure concepts to demonstrate its applicability in both highly equipped canopy research facilities and forests without permanent canopy infrastructure. One installation relied on a canopy crane, grid power and fully automated irrigation, whereas the second was built by certified tree climbers and was equipped with an autonomous solar-powered, battery-operated irrigation system. At both sites, environmental sensor networks monitor the experiment. 4. TreeTOP successfully reproduced characteristic canopy microclimatic gradients, including increasing light availability, daytime air temperatures and thermal extremes with canopy height. Despite differing infrastructures, both implementations generated comparable microclimatic patterns, demonstrating that standardized canopy experiments are feasible in forests with or without permanent canopy access. By opening canopy space for manipulative experiments, TreeTOP provides a transferable framework for investigating plant performance, phenology, species interactions and microbiome assembly under realistic forest conditions.

Citizen science observations are abundant, but conservation requires turning uneven records into reliable predictions and directing new surveys to where information is missing. We developed a biodiversity platform for Japan that is updated monthly and integrates 2.32 million records to predict 8,297 species across seven taxonomic groups. Shared representation models outperformed species-specific models in four groups and extended predictions to species with few records. Five independent datasets, including structured monitoring, environmental DNA and complete forest inventories, confirmed that the models ranked observed species and occupied sites above alternatives, with median AUCs of 0.724 to 0.894 across sites and 0.650 to 0.841 across species. For any user-selected area, the platform returns candidate species, distribution predictions, a biodiversity map corrected for uneven observation effort, a conservation priority map for native species and a map recommending where to survey next. This map highlights places where species with few records are predicted to occur despite limited sampling. Independent observations showed that areas ranked highly by this predicted potential contained many such species, indicating that model predictions can help direct surveys toward knowledge gaps. New observations are incorporated into monthly updates, creating a national feedback system connecting citizen science, local conservation decisions and future surveys.

Heritable bacterial symbionts are pervasive in terrestrial arthropods, often imposing reproductive manipulations to promote their own spread within host populations. Co-infections are common, potentially allowing symbiont co-infectors to hitchhike through a host population. However, adverse thermal conditions can disrupt these communities, particularly when co-infectors vary in their thermal sensitivity. We used a multi-generation experiment to test whether warm (29 {degrees}C) conditions disrupted spread of heritable symbionts through uninfected populations of the spider, Mermessus fradeorum. We tested two common infection combinations: a single infection with a cytoplasmic incompatibility (CI) inducing Rickettsiella or a feminizing co-infection that included a feminizing Wolbachia, the same Rickettsiella, and up to three apparent hitchhikers (two additional Wolbachia strains and Tisiphia). We initiated replicate populations with 1/3 of one infection type and 2/3 uninfected spiders, evaluating population infection rate over 5 spider generations under different temperature regimes. Under cool (21{degrees}C) conditions, Wolbachia feminization drove co-infection to 88% and Rickettsiella CI drove single infection to 83% of host populations. Vertical transmission for all symbionts was high (97-99%) and hitchhiking symbionts also spread effectively. Under warm conditions, feminization and CI efficacy were reduced, and symbionts suffered variably reduced vertical transmission. Warm conditions ultimately destroyed the co-infecting symbiont consortium and impeded symbiont spread. On its own, though, Rickettsiella was still able to increase, despite reduced strength of CI. We hypothesize that contrasting tensions between feminizing spread of the symbiont consortium versus environmentally driven loss of function and transmission may explain observed patterns of mixed infections in field populations of this spider.

Plant cold acclimation emerges from coordinated adjustments in photosynthesis, primary metabolism, and intracellular carbon allocation. Yet, the regulatory role of subcellular metabolite compartmentation in natural variation of cold acclimation remains insufficiently understood. Here, we investigated four Arabidopsis thaliana accessions grown either individually or in bulk to determine how growth configuration and genotype shape the metabolism of sugars and organic acids during cold exposure. Using non-aqueous fractionation, we quantified plastidial, cytosolic, and vacuolar sugar pools alongside whole-cell carbohydrates, organic acids, enzyme activities, photosynthetic parameters, and stress markers. A neural-network classifier revealed that subcellular sugar distribution together with sugar amounts and organic acids provided the strongest discriminatory power among accessions, surpassing photosynthetic traits and enzyme activities. Our findings demonstrate that natural variation in cold acclimation is strongly determined by genotype-specific subcellular metabolite architectures, and that the cultivation strategy modulates these intracellular signatures. We conclude that subcellular compartmentation of metabolites represents a cellular control layer for natural variation of cold acclimation and resilience in Arabidopsis thaliana.

Standard methods for ancestral sequence reconstruction (ASR) rely on substitution models for the residues in a biological sequence and assume independent evolution across these sites, ignoring the epistatic interactions that shape molecular evolution. In contrast, deep learning models like variational autoencoders (VAEs) can learn low-dimensional representations ("embeddings") of sequences in a protein family that may implicitly handle these dependencies, raising the possibility of performing more accurate ASR by interpolating between extant sequence embeddings within the VAE's latent space. In this study, we test this hypothesis by developing and evaluating a VAE-based ASR pipeline. Benchmarking this approach against established likelihood-based and parsimony methods using various simulations of protein evolution, including scenarios with and without epistasis, we find that the VAE-based approach is consistently and significantly outperformed by standard methods, even in epistatic regimes where it was hypothesized to have an advantage. We further show that this failure is not due to a lack of phylogenetic structure in the latent space, which does contain evolutionary signal. Rather, the primary limitation is the information loss inherent to the autoencoding process: the VAE's decoder cannot generate sequences with sufficient fidelity for the precise demands of ASR.

Infections with bacteria which have evolved multi-drug resistance (MDR) cause millions of deaths worldwide. Large-scale efforts are gathering genotypic and phenotypic data on MDR bacteria, but methods for learning the structure, diversity, and predictors of evolutionary pathways to MDR have yet to take full advantage of these data. Here, we use evolutionary accumulation modelling (EvAM), an emerging class of machine learning methods with roots in cancer progression, to infer these evolutionary pathways across ESKAPEE pathogens (seven bacterial species that dominate health burdens), using a database of over 635k genotyped phenotypic observations from around the world. We identify global patterns in MDR evolutionary pathways, remarkably shared across multiple ESKAPEE species. Species-specific deviations from these stereotypical pathways are connected with geographical and demographic covariates, facilitating predictions of future MDR evolution. We verify these predictions with several hundred new phenotypes from ESKAPEE samples spanning decades of clinical infections in sub-Saharan Africa, demonstrating the capacity to forecast future MDR evolution from these inferred shared pathways.

Marine diatoms contribute to about 20% of global primary production. We present the first basin-scale, multiyear assessment of diatom communities in the North Pacific, combining taxonomically high-resolution RuBisCO large subunit gene (rbcL) metabarcoding with concurrent environmental measurements. Using a nine-year time series of daily samples resolved at the species level via ~500 bp rbcL fragments, we performed multivariate analyses across biogeographic provinces, identifying significant correlations between community structure and environmental drivers such as temperature and macronutrient availability. We report the prevalence of a previously overlooked centric diatom species in the North Pacific, Eunotogramma lunatum, which appears to be near-dominant even in subarctic high-nitrate, low-chlorophyll waters where pennate diatoms are typically favored. These results demonstrate the power of rbcL for large-scale ocean monitoring and provide a critical baseline for future studies of diatom population dynamics, climate change impacts, and ecosystem resilience in a key marine region.

Phage tail-like elements (PTEs) -- tailocins, bacterial type VI secretion systems (T6SS), and extracellular contractile injection systems (eCIS) -- are contractile nanomachines that bacteria use to kill their neighbors and compete within their micro-ecosystems. PTEs help shape microbial community composition. Most PTE detection tools only detect a single PTE class. Moreover, most tailocin detection methods are largely restricted to Pseudomonas, leaving a key part of tailocin diversity uncharacterized. In this work, we present PhageTAILor (https://github.com/hjcho-bio/PhageTAILor), an integrative and fully automated pipeline that detects and classifies prophages and 3 PTE classes from bacterial genomes. PhageTAILor combines a 6-detector homology-based candidate search (geNomad, tail-gene, PHROGs-tail, SecReT6, eCIStem, and a divergence-tolerant tail-HMM detector) with a LightGBM classifier comprising 1 multiclass and 3 binary heads, trained on 6,501 bacterial genomes carrying 13,082 prophages and PTEs. A phylogeny-free feature matrix used in our model keeps predictions reproducible between model construction and user inference. PhageTAILor performs strongly at the genome level and generalizes beyond its Pseudomonas-rich training set. On a 76-strain cross-clade benchmark, PhageTAILor detected tailocins at F1 = 0.955. Furthermore, it identified 12 of 13 experimentally validated tailocins spanning five genera versus 2 of 13 for a Pseudomonas-restricted tool TattleTail. PhageTAILor also demonstrated sensitivity equivalent to viral detection tool geNomad while avoiding its higher false-positive rate. Applied to 7,925 plant- and soil-associated bacterial isolates, PhageTAILor showed that prophages in the phyllosphere and tailocins in plant-associated bacteria, whereas eCIS are enriched in soil. PhageTAILor is distributed as an open-source, modular pipeline with a command-line interface.

Across different modalities, signals play a core role in attracting mates and influencing mating success. In several non-human primate species, females produce calls during mating that are thought to promote male competition over receptive females. The extent to which social system characteristics modulate the function of copulation calls remains less clear. We studied copulation calls in wild Guinea baboons (Papio papio), who live in a multilevel society structured around units in which females associate and mate almost exclusively with a single male. We hypothesised that females use copulation calls as an indirect form of mate competition, with competition increasing in larger units. In addition, we hypothesised that females are more likely to mate again after calling. We analysed 6116 copulations between 2014 and 2025, involving 99 reproductively active females and 78 subadult and adult males. Females produced copulation calls in 72.7% of copulations, with large inter-individual variation. Neither unit size nor its interaction with the female's swelling size or the presence of simultaneously receptive females affected the probability of calling. A survival analysis with a subset of the data (2353 copulations) revealed no effect of calling on the latency to the next mating. Our results render the hypothesis that female Guinea baboons use calls in indirect mate competition unlikely. Yet, the probability of calling varied with sexual swelling size, suggesting that calls signal female fertility. Possibly, Guinea baboon copulation calls represent an evolutionary remnant, no longer under selective pressure, and can be considered index signals of female fertility.

All living organisms rely on the movement of ions across cell membranes as the fundamental physical basis of their internal energy and signaling, and plants are no exception. Plants perceive, integrate, and respond to environmental stimuli through electrical signals, classified as action, variation, and system potentials, that are coupled with calcium waves, reactive oxygen species, and hydraulic and hormonal changes to coordinate whole-organism responses despite the absence of a nervous system. Yet most studies characterize these signals using a single feature, such as amplitude or spike duration, in a single tissue, an approach that cannot establish how such signals correspond to the underlying ionic activity, mobility, and structural complexity of the signaling environment, or how this correspondence varies across organs. Here, we correlate plant bioelectrical signals with potential ionic energy flow using a multi-domain framework, combining discrete spike events, continuous waveform properties, spectral composition, and signal complexity applied to leaf, stem, and root recordings from tomato (Solanum lycopersicum) exposed to different stimulus. Electrical activity with increased stimulus strength, likely reflecting increased ionic flow, with the root showing the largest response. This suggests plant electrical signaling works as a distributed, ion-based information system, useful for stress monitoring and bio-inspired sensor design.

Wildlife re-identification (re-ID) is a widely used and powerful tool with diverse applications in animal ecology and conservation. Current automated methods typically operate on single images of a single body part of the animal. However, a single encounter may contain multiple images capturing different body regions, each providing complementary individual-specific information. In contrast to automated approaches, researchers often manually select the most suitable images and regions for identification based on factors like visibility, occlusion and image quality. This creates a mismatch between automated methods and field practice, limiting the practical adoption of current automated re-ID pipelines. Here, we address this by introducing an encounter-based, multi-body-part re-ID framework, using sea turtles as a model taxon. Our framework combines three elements: (1) An orientation-aware deep learning model, TurtleDetector, that in addition to the full bodies, it also automatically segments key body regions, i.e. heads, front and hind flippers, from images within an encounter; (2) a hybrid body-part-specific retrieval method, that sequentially combines a fast global-feature model (MiewID or DINOv3) with a more accurate but costlier local-feature model (ALIKED with LightGlue); and (3) a merged identity-prediction strategy that selects the highest calibrated similarity score across all available body parts and images of an encounter. We evaluate the framework on three long-term re-ID datasets spanning three species, loggerheads, greens, and hawksbill turtles, under an evaluation protocol that mirrors real-world, time-aware re-ID workflows. Across datasets, combining multiple body regions consistently improved identification performance over the best-performing single body region, resulting to an increase of 4-6% in top-1 accuracy. Interestingly, body regions traditionally underused in sea turtle re-ID, such as the hind flippers and carapaces, provided complementary identifying information that improved encounter-level re-ID when integrated through the hybrid retrieval method. Our findings demonstrate that automated wildlife re-ID can benefit from moving beyond single-image, single-body-part identification towards encounter-level integration of all available visual evidence. Our work further suggests that, where feasible, field photo-acquisition protocols should aim to capture multiple informative views of an individual during each encounter. Importantly, many species and taxa, including elephants, primates, cetaceans, and other large vertebrates, possess such individual-specific features across multiple body regions, highlighting the broad potential applicability of our framework.

DNA sequencing unified microbial genotyping into a single, comprehensive readout, yet phenotyping remains a slow and fragmented endeavor. Here, we introduce Microbial Phenotyping Using Low-magnification Label-free Imaging (PULLI), a computer vision platform that extracts microcolony and population-level phenotypes from brightfield timelapses of liquid culture growth. Using PULLI, we screened a genome-scale Vibrio cholerae mutant library, recording more than 200,000 images, which revealed that core bacterial pathways shape community architecture. Functionally related mutants converge in appearance, allowing us to resolve processes as distinct as biofilm formation, motility, central metabolism, cofactor biosynthesis, and envelope composition using a single approach. We further show PULLI can be used to determine a drug target, characterize other pathogens, and classify bacterial species. Our results show that bacterial multicellular development is an interpretable signature of genotype-phenotype relationships, which can be captured from simple brightfield timelapses. We release the PULLI pipeline and an interactive atlas of community forms.